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Secondary electron emission due to Auger de-excitation of 
metastable nitrogen molecules at metallic surfaces 
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With an eye on plasma walls we investigate, within an effective model for the two active 
electrons involved in the process, secondary electron emission due to Auger de-excitation of 
metastable nitrogen N2(^Ti^) molecules at metallic surfaces. Modelling bound and unbound 
molecular states by a LCAO approach and a two-center Coulomb wave, respectively, and the 
metallic states by the eigenfunctions of a step potential we employ Keldysh Green's functions 
to calculate the secondary electron emission coefficient 7e retaining for the Auger self-energy 
the non-locality in time and the dependence on the single particle quantum numbers. For 
the particular case of a tungsten surface we find good agreement with experimental data 
indicating that the relevant Auger physics is well captured by our easy-to-use model. 



1. Motivation 

Secondary electron emission by surface scatter- 
ing of metastable molecules is an important pro- 
cess in molecular low-temperature gas discharges. 
It is one of the main wall-based electron gen- 
eration channels controlling, together with wall 
recombination and various volume-based charge 
production and destruction channels, the overall 
charge balance in the discharge. The efficiency of 
electron release can be encapsulated in the sec- 
ondary electron emission coefficient 7e, which de- 
notes the probability for releasing an electron by a 
single metastable de-excitation at the wall. How- 
ever, for most wall materials and projectiles 7e 
is unknown. A flexible, easy-to-use microscopic 
approach for its calculation is thus needed. 

Stracke et al. [1] experimentally investigated 
the de-excitation of metastable nitrogen N2{^Tj^) 
molecules on a tungsten surface and proposed two 
reaction channels. Firstly, Auger de-excitation 
(also referred to as Penning de-excitation when 
the emitted electron comes from the molecule), 
A^2(^S+) + em ^ ^^2(^S+) + ef, and, secondly, 
the formation of a Ni^C^Ug) shape resonance with 
subsequent auto-detachment, N2{^'S^) + Cm ^ 
N2(^^g) — > N2(}Ti'^) + e/, where em and e/ de- 
note an electron inside the metal and a free elec- 
tron, respectively. Stracke et al. [1] concluded 
that the latter process should be more efficient, as 
it is a combination of two single-electron charge- 
transfer transitions, whereas the Auger reaction 
represents a less probable two-electron transition. 
Using thermal molecules they measured the en- 
ergy spectrum of the released electron and esti- 



mated 7l^''|totai to be about 10"^ — 10~^. This 
value does however not discriminate between the 
two channels. Indeed, in the energy-resolved data, 
Stracke et al. [I] also find a weak signal due to 
Auger de-excitation giving rise to 7!^'' | Auger ~ 
10~^ — 10~^. In view of the electronic structure of 
tungsten this is not surprising because large en- 
ergy shifts and broadening due to image interac- 
tions are required to bring the molecular orbital 
hosting the hole in the electronic configuration 
of N2(^Ti'^) in resonance with conduction band 
states. Charge transfer may thus be less efficient 
and Auger de-excitation an eye-to-eye competitor. 

To establish our approach and because of 
the availability of experimental data we focus 
in the present work on Auger de-excitation of 
iV2(^S^) at an uncharged tungsten surface with 
simultaneous release of an electron. We construct 
an effective model that concentrates on the most 
important degrees of freedom and enables us to 
describe the system by a few parameters which 
are accessible through experiments or theoretical 
calculations. The primary goal will be to calcu- 
late 7e using Keldysh Green's functions along the 
lines layed out by Makoshi [2] in his investigation 
of de-excitation of metastable atoms. 

2. Model 

Focussing on the essentials of the process, we 
assume the metallic surface to be planar, ideal, 
and to stretch over the entire half space z < 0. 
Furthermore, we consider only the dominant 
metastable state N2(^Ti'^) and employ the tra- 
jectory approximation, that is, we decouple the 
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translational motion of the molecule from the dy- 
namics of the system and externally supply its 
trajectory. Finally, the molecule is assumed to 
impact the surface under normal incidence with 
constant velocity v and constant angle ip of its 
axis to the surface. Assuming the molecule to 
start moving at to = -co and to hit the sur- 
face at t = the trajectory of its center of mass 
is R{t) = {y\t\ + zq) ez 1 where zq is the turn- 
ing point of the molecule in the surface potential 
Vs{z). We consider only indirect de-excitation 
(Penning process) in which the electron is emit- 
ted from the molecule (solid lines in Fig. [T|). Di- 
rect de-excitation (dashed lines in Fig. [T|) can be 
shown to be negligible because the orthogonality 
of the bound molecular wave functions makes the 
corresponding matrix element much smaller than 
the one for Penning de-excitation. 

To construct the Hamiltonian we combine three 
different kinds of single-electron states to a single- 
electron basis: the single-electron states of the 
conduction band of the solid surface \k), which 
we approximate by the states corresponding to 
an electron trapped by a step potential of depth 
$(7) the free single-electron states \q} associated 
with the molecule's continuum for which we use a 
two-center Coulomb wave, j3j and effective single- 
electron states |Om/lm) (m = ±1 is the magnetic 
quantum number) for the bound states of the 
molecule, which we approximate by a degenerate 
two-level system keeping, within the LCAO rep- 
resentation of the nitrogen molecule, only the 277^ 
and the 27rg molecular orbitals (MOs) which are 
the two MOs whose occupancies change during 
the de-excitation process. The LCAO MOs are 
constructed by hydrogen-like wave functions with 
effective nucleus charges to mimic the Roothaan- 
Hartree-Fock wave functions of atomic nitrogen. 
Note, the Penning process does not involve any 
spin flip, we thus ignore the spin. 

The description of the electronic structure of 
the molecule-surface system requires to align 
the single-electron states against each other and 
against the vacuum level by use of the metal's 
work function the metal's conduction band 
depth $(7, the molecule's ionization energy Ae^, 
and the excitation energy of the molecule Asg- 
The metal states are of course occupied up to the 
Fermi level ep- Our model is thus characterized 
by a few energy parameters, an effective charge, 
and a bond length (which enters the molecular 




Fig. 1: Illustration of the effective model showing Pen- 
ning de-excitation (solid lines) and its exchange pro- 
cess (dashed lines). Here e^*" and e'^'^^ are the clas- 

k y 

sical energy cut-offs obtained from the energy balance 
£l -\- = £q -\- £^ that holds in the adiabatic limit. 

wave functions). 

Due to the symmetries of the molecular ground 
and excited state, only transitions with Am = 
are involved. They are driven by the Coulomb 
interaction between the excited electron in the 27rg 
MO and an electron in the Fermi sea of the metal. 
The three electrons in the 27r„ MO act only as 
spectators and can thus be neglected. Assuming 
moreover the Fermi surface of the metal to be 
rigid, the de-excitation of A^2(^5]+) is a two-body 
scattering process, whose Hamiltonian, written in 
the single-electron basis just described, is given by 

k k k ^ gig 

k 9 

+ ^ ^0 ^Om ^Om + ^ ^1 '^Im ^Im (1) 
m m 

k,q,m 

The non-equilibrium character of the system, and 
hence the need to employ Keldysh Green's func- 
tions [2], arises from the time-dependent Penning 
matrix element, 

K!;:r(*) = jdrjdr' %^{r^(t))^.{r) 

X Vc{\r-r'\)^l{r'^{t))^,^X^'^{t)) , 

(2) 
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where Vc represents the Coulomb potential and 
the subscript ip denotes the associated vector as 
seen from the molecule's reference frame. 

Image interactions are not considered except 
for the emitted electron, which always feels 
the image potential Vi{z) and can only escape 
from the surface when its perpendicular energy 
is sufficiently high. To take this into account 
we multiplied in the calculation of je and the 
spectrum of the emitted electron the Penning 
matrix element by the surface transmission 
function Q{eq^ - Vi{zR{t))). 

3. Quantum kinetics 

The secondary electron emission coefficient can 
be obtained from the time-dependence of the oc- 
cupancy of the free states \q). To calculate this 
quantity from dH) we employ Keldysh Green's 
functions along the lines pioneered by Makoshi. [2j 
In contrast to him, however, we do not work with 
a phenomenological Auger interaction and do not 
employ the wide-band approximation for the free 
states. We are also not restricted to the lowest or- 
der expressions for the occupancies (given below 
by (HI), dS]), and but can in principle calculate 
higher order corrections to these expressions. 

The central quantity of our approach (in detail 
described in Ref. 0]) is a function Am(ti,t2) 



A, 



q,k 



(3) 



which emerges from the self-energy terms of the 
iterated Dyson equation. To lowest order in 
the occupation of the excited molecular level is 



(4) 



with fi being the initial magnetic quantum num- 
ber of the excited electron. For the spectrum of 
emitted electrons we obtain to lowest order in A^ 

ng-(t) = fdh f dt2 A^,,-(ti,t2) , (5) 

where A^ g- is defined by Eq. ([3]) without the sum 
over q. 

The secondary electron emission coefficient 7e 
can be calculated from Eq. ^ by taking t = oo 
and summing over all possible q, which to lowest 
order in A^, yields 



7, 



(0) 



poo poo 

/ dti / dt2 A^.(ti,t2) , 
Jto Jto 



(6) 



where fi is again the initial magnetic quantum 
number. If we allow for the next order terms with 
respect to A^ in ([5]) we find 



7, 



(1) 



1 



(7) 



However, the zeroth and first order results ([6]) and 
([7]) differ significantly only for very small kinetic 
energies of the molecule. [3] 

To simplify the numerical analysis of ©-([T]) 
we utilize the localized nature of the molecular 
wave functions, which together with the fact, that 
the molecule's turning point zq lies far outside 
the surface (zq > 4.35 as) allows us to neglect the 
oscillating part of the metal wave functions inside 
the solid and to separate the time dependence 
of the Penning matrix element. [1] The time 
integrals of Am{ti,t2) can then be calculated 
analytically. The remaining wave vector depen- 
dent part of the matrix element is calculated 
on a discrete grid within the {k, ^)-space and 
afterwards the wave vector integrals within A^ 
are calculated using sexa-linear interpolation to 
obtain the inter-grid values of V'-^"!, 

4. Results 

We now present selected results for the particu- 
lar case of a tungsten surface. The kinetic energy 
of the molecule was fixed to 50meV. The associ- 
ated turning point zq amounts to approximately 
4.4:2 as- We did not find any influence of the ini- 
tial magnetic quantum number = ±1. For con- 
venience we thus omit any m subscripts in the 
following. In addition we restrict the molecule's 
orientation to the two fundamentally distinct sit- 
uations if = (axis parallel to the surface) and 
93 = ^ (axis perpendicular to the surface). 

Figure [2] shows the time evolution of the occu- 
pancies of the excited molecular level and the free 
electron states. Obviously, the occupancy of the 
excited molecular level (upper panel) changes sig- 
nificantly only in the range \t\ < 5, which, taking 
the turning point into account, equals maximum 
distances of the molecule's center of mass from the 
surface of roughly 7aB- Within our model, the 
process is equally effective in the incoming and 
outgoing branch of the trajectory. The time evo- 
lution of the occupancy of the free states (lower 
panel) is however distinctively different. It shows 
a plateau around t = 0, that is, a stagnation of 
the probability to escape from the solid. This is a 
consequence of the image potential which almost 
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Fig. 2: Time evolution of the occupancies of the ex- 
cited molecular level and the free electron states for 
parallel (solid lines) and perpendicular (dashed lines) 
molecule orientation. 




Fig. 3: Spectrum of the emitted electron for parallel 
(solid lines) and perpendicular (dashed lines) molecule 
orientation. 



completely traps an electron emitted at low sur- 
face distances where its perpendicular energy is 
too small to overcome the image barrier encoded 
in the surface transmission function. The number 
of emitted electrons is larger in the perpendicular 
case as compared to the parallel case. Numerical 
tests showed that for the parallel case the electron 
is primarily emitted with a very small perpendic- 
ular energy whereas in the perpendicular case the 
perpendicular energy of the emitted electron is 
distributed more equally. Thus, in the perpendic- 
ular orientation the electron has a higher proba- 
bility to breach through the surface barrier origi- 
nating from the image potential. The occupancy 
of the free states for t 00, that is, the secondary 
electron emission coefficient amounts to 2 • 10"'^ 
in the perpendicular orientation and 2.4 • 10~^ in 
the parallel orientation which is quite close to the 
experimental estimate given by Stracke and coau- 
thors, I Auger ~ 10-^ - 10-3. [T] 

Figure [3] finally shows the spectrum of the 
emitted electron at t = 00. The graphs for the 
two different orientations start at the origin and 
monotonously increase until a cut-off energy 
is reached. The latter resembles the classical 
cut-off energy e'^'^^ (see Fig. [1]), implying that 
energy conservation is restored at the end of the 
collision, as it should be. The low energy part of 
the spectrum is cut off due to the surface trans- 
mission function which allows electrons to escape 
from the surface only when their perpendicular 
energy is large enough. The spectrum for the 
perpendicular case takes on larger values, for 
reasons already explained. 



5. Conclusions 

We investigated within the Keldysh formalism 
the release of secondary electrons due to Auger 
de-excitation of N2{^T,^) at metallic surfaces 
using an effective model for the two active 
electrons involved in the process. For tungsten 
we obtained for this particular electron emission 
channel 7e ~ 2 • IQ-^ (perpendicular orientation) 
and 7e ~ 2.4 • 10"^ (parallel orientation) which 
agrees well with experimental estimates. The 
effective model depends only on a few, easily 
obtainable energies - a real advantage for appli- 
cations to low-temperature gas discharges, where 
a great variety of different kinds of molecules and 
wall materials occur. With some modifications 
the model can be applied to dielectric surfaces 
as well. We employed the model only for the 
investigation of Auger de-excitation. But pro- 
vided the molecule's image interactions and the 
on-molecule Coulomb interaction between the 
two active electrons are taken into account, it can 
be also used to study charge-transfer processes. 
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